A1022 


•A 


Qualified  requestors  may  obtain  additional  copies  from  the 
Defense  Technical  Information  Center.  All  others  should 
apply  to  the  National  Technical  Information  Service. 


SECURITY  CLASSIFICATION  of  this  *>AOE  (Wht t  Timtm  Entered) 


REPORT  DOCUMENTATION  PAGE 


RF>  r>  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


.  REPORT  NUMDER 


AFGL-TR-80-0247 


4.  TITLE  (m«t  SublllU) 


Satellite  Data  Processing  Systems  and 
Accelerometer  Density  Data  Base 


7.  AUTHORS) 

Robert  W.  Fioretti 
Andrew  J.  Mazzella,  Jr. 

Rebecca  J.  Leong,  Lawrence  D.  Cox 


19.  performing  organization  name  and  aooress 


CATALOG  NUMBER 


5.  Tyre  OF  KEPC/AT  A  PERIOD  COVERED 

1  Sept.  78  -  31  Mar.  80 


6  PERFORMING  ORG.  REPORT  NUMBER 

FINAL  REPORT 


0.  CONTRACT  OR  GRANT  NUMBER*'*) 

F19628-78-C-0226 


10  PP'  C.PAM  PL  fMENT.  PROJECT.  TASK 
4»PI  a  A  WORK  UNIT  NUMBERS 

A  *  A  62101F 

391  Totten  Pond  Road  ccon  g-jbc 

_ Wal  ■Fham  f  MA  Q2I5.4 - -  - 

II.  CONTROLLING  OFFICE  NAME  AND  AODRESS  '*•  REPORT  DATE 

Air  Force  Geophysics  Laboratory  30  April  1980 

Hanscom  AFB,  MA  01731  is.  njmbeVoV  pages 

Monitor/Frank  Marcos/LKB _ 31 _ 

TT  MONITORING  AGENCY  NAME  a  ADORESS(l/  different  from  Controlling  Office)  15.  SECURITY  CLASS,  (of  thle  report) 

UNCLASSIFIED 


30  April  1980 

13.  N'JMBEpTTf  PAGES 


15*.  OECL  ASSf  FI  CATION/ OOWNGR  AOING 
SCHEDULE 


I  l«.  DISTRIBUTION  STATEMENT  (of  this  Report) 


Approved  for  public  release;  distribution  unlimited. 


I  17.  DISTRIBUTION  STATEMENT  (of  the  mbmtemct  entered  In  Block  20,  U  different  from  Report) 


10.  SUPPLEMENTARY  NOTES 


19.  KEY  WOROS  (Continue  on  revetee  eld#  II  neceeemry  end  Identity  by  block  number) 

Data  Processing,  Accelerometer,  Density,  Data  Base, 
Regression  Modelling. 


I  20.  ; ABSTRACT  (Continue  on  I 


Hide  If  neceeeery  mod  Identity  by  block  number) 


This  report  describes  the  finalization  of  an  extensive  neutral 
atmospheric  density  data  base  taken  from  four  satellite  accel¬ 
erometer  systems;  development  of  an  empirical  model  of  these 
data,  circular  orbit  processing,  AE-C  re-entry  data  and  contrac¬ 
tual  functions  in  support  of  the  EUVS  experiment. 


DD  ,:sr»  1473  EDITION  OF  1  NOV  SS  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  r**»n  Jimtm  Entfdl 


SECURITY  CLASSIFICATION  OF  THIS  PAQEflFhan  Data  Bntarad) 


SECURITY  CLASSIFICATION  OF  THIS  PAGEflWiwi  Data  Bntarad) 


FOREWORD 


The  efforts  described  herein  were  performed  under 
contract  to  the  Atmospheric  Structure  Branch  (LKB) , 
Aeronomy  Division  of  the  Air  Force  Geophysics  Labora¬ 
tory  (AFGL) ,  Hanscom  Air  Force  Base,  Massachusetts. 


'  viVut :  -::i  /  | 

Avaiia'.-ililv  j 

an.:/ or  j 

Diet  j  Special  j 


Table  of  Contents 


Foreword .  iii 

Table  of  Contents .  v 

List  of  Illustrations .  vi 

1.  Introduction .  1 

1.1  AE/S3  Data  Base .  1 

1.2  New  AE-C  Orbits .  3 

1.3  Data  Base  Properties .  3 

2.  Model  Development .  6 

2.1  J71,  J77  and  MSIS  Comparison..  6 

2.2  Empirical  Model .  6 

2.3  Results . 12 

3.  Other  Contract  Functions . 19 

3.1  EUVS  Experiment  Support . 19 

3.2  AE-C  Re-Entry  Data . 21 

3.3  Circular  Orbit  Processing . 21 

4.  Conclusions . 22 

5.  Acknowledgment . 22 

6.  References . 2  3 


v 


priEC*I)UO 


List  of  Illustrations 


1  Distribution  of  Density  Data .  2 

2  Data  Distribution  as  a  Function  of 

Geomagnetic  Latitude .  4 

3  Data  Distribution  as  a  Function  of 

Local  Time .  4 

4  Data  Distribution  as  a  Function  of 

Geomagnetic  Activity .  4 

5  Data  Distribution  as  a  Function  of 

Solar  Flux .  5 

6  Data  Distribution  as  a  Function  of 

Altitude .  5 

7  Data  Base  at  180  km:  Solar  Flux, 

Density,  Ratio  to  J71  Model 

and  Kp  as  a  Function  of  Time..  7 

8  Frequency  Distribution  of  Percent 

Differences  of  all  Data 

from  J71  Model .  8 

9  Frequency  Distribution  of  Percent 

Differences  of  all  Data 

from  MS  IS  Model .  8 

10  Frequency  Distribution  of  Percent 

Differences  of  all  Data 
from  J77  Model .  8 


] 1  Frequency  Distribution  of  Percent 

Differences  of  200  km  Data 
from  J71  Model . 

12  Frequency  Distribution  of  Percent 

Differences  of  200  km  Data 


from  MSIS  Model .  9 

13  Frequency  Distribution  of  Percent 

Differences  of  200  km  Data 

from  J77  Model .  9 

14  Amplitude  and  Phases  of  Annual  and 

Semi-Annual  Variations  of 

the  Data  Base..... .  13 


List  of  Illustrations  (Cont.) 


Figure  page 


15  Amplitude  and  Phases  of  Diurnal  and 

Terdiurnal  Variations  of 

the  Data  Base .  14 

16  Semi-Annual  Density  Variation  for 

each  Model  at  200  km  for  0  Degrees 
and  60  Degrees  Latitude .  15 

17  Geomagnetic  Activity  Variations  of 

Density  for  each  Model  at  200  km 

for  0  Degrees  and  60  Degrees 

Latitude .  16 

18  Solar  Flux  Variations  of  Density  for 

each  Model  at  200  km  for  0  Degrees 
and  60  Degrees  Latitude .  17 

19  Model  Density  Fit  Results  as  a  Function  of 

Altitude  for  Days  80  and  172 .  18 

20  Density  Contours  as  a  Function  of  Latitude 

and  Day  for  Ap  =  0 .  20 

21  Density  Contours  as  a  Function  of  Latitude 

and  Day  for  Ap  =  80 .  20 


vii 


> 


1.  INTRODUCTION 


The  efforts  described  herein  are  part  of  a  project 
to  create,  maintain  and  utilize  an  extensive  neutral 
atmospheric  density  data  base  qenerated  from  four  AFGL 
satellite  accelerometer  experiments  -  three  on  the  NASA 
Atmosphere  Explorer  (AE)  series,  AE-C,  AE-D,  AE-E,  and 
one  on  the  Air  Force  S3-1  satellite.  The  density  data 
were  taken  over  a  three  year  period  from  December  1973 
to  November  1976.  The  accelerometer  data  processing  sys¬ 
tem,  the  analytical  techniques  used  to  extract  atmospheric 
density  from  the  total  accelerometer  output  signal,  and 
the  creation  of  a  density  data  base  have  been  described 
previously  (1,2,3). 

This  report  describes  the  finalization  of  an  extensive 
neutral  atmospheric  density  data  base  taken  from  the  four 
satellite  accelerometer  systems,  the  development  of  an  em¬ 
pirical  model  of  these  data,  and  other  aspects  of  data 
analysis  support  for  AFGL. 

Studies  by  AFGL  scientists  utilizing  these  results  are 
described  elsewhere  (e.g..  Ref.  (4)). 

1.1  AE/S3  Data  Base 

The  AE/S3  density  data  base  includes  data  from  the  ellip¬ 
tical  orbit  portions  of  the  respective  missions  of  the  AE-C, 
-D,  -E  and  S3-1  satellites.  Table  1  provides  orbital  char¬ 
acteristics  and  dates  of  elliptical  orbit  data  acquisition 
for  each  satellite.  The  altitudes  at  which  data  were  obtained 
are  shown  in  Figure  1. 


1 


Table  1.  Satellite  Orbital  Characteristics  and  Data  Acquisition  Periods 


Satellite 

Launch 

End  Elliptical 
Data  Acquisition 

Inclination 

AE-C 

'  Dec  73 

Nov  74 

68° 

S3-1 

Oct  74 

May  75 

97° 

AE-D 

Oct  75 

Jan  76 

90° 

AE-E 

Nov  75 

Nov  76 

20° 

Perigee  Apogee 

156  km  4000  km 
160  km  4000  km 

156  km  3800  km 

157  km  3000  km 


eseei 


50  100  150  200  Z50  300  350 1  50  KX)  150  200  250  300  350 1 
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Figure  1.  Distribution  of  Density  Data 


Density  values  were  calculated  at  5-km  intervals 
and  stored  in  the  data  base.  Associated  with  each  point 
is  the  corresponding  satellite  enhemeris  data,  altitude 
information,  solar  flux,  Kp/Ap  values,  J71,  J77  and 
MS  IS  atmospheric  models.  References  (3)  and  (4)  describe 
this  data  base  in  detail- 
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1. 2  New  AE-C  Orbits 

More  than  one  thousand  (1000)  orbits  of  AE-C 
elliptical  orbit  data  not  previously  processed  became 
available  and  were  processed  under  the  present  contract. 
These  data  were  processed  using  the  MESA  data  reduction 
system  (previously  described  in  Reference  (1)),  and 
calculated  density  values  were  stored  in  MESA  GU  and 
UA  files. 

The  new  AE-C  density  data  base  was  then  developed 
from  the  new  UA  density  data  using  the  techniques  des¬ 
cribed  in  Reference  (3).  This  data  base  was  generated 
on  the  Sigma-9  computer  and  was  written  to  magnetic  tapes. 
These  tapes  were  shipped  to  AFGL  for  insertion  into  the 
AE/S3  data  base  on  the  AFGL  CDC  6600  computers. 

To  accomplish  this,  programs  were  developed  to  insert 
(i.e.,  merge)  the  new  data  into  the  existing  AE-C  data 
base  for  all  altitudes,  to  quality  check  the  resultant 
merged  data  base  for  duplicate  data,  and  finally  to  re¬ 
merge  this  data  base  with  the  remainder  of  the  AE/S3  data 
at  all  altitudes.  This  new  composite  data  base  makes  up 
the  present  AE/S3  density  data  base. 

1. 3  Data  Base  Properties 

Figures  2-6  are  histograms  of  the  data  distribution 
as  a  function  of  geomagnetic  latitude,  local  time,  Kp 
index,  solar  flux  and  satellite  altitude.  Figure  2  shows 
the  number  of  data  points  obtained  in  each  20  degree  geo¬ 
magnetic  latitude  band.  Figure  3  shows  the  distribution 
of  data  points  obtained  for  each  hour  of  local  time.  The 
distribution  of  data  with  respect  to  geomagnetic  activity 
is  given  in  Figure  4.  This  gives  the  number  of  measure¬ 
ments  during  which  the  three-hourly  Kp  index  (with  6  hr. 
lag)  fell  within  a  soecific  one  unit  interval.  Tn  Figure  5 
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Figure  2.  Data  Distribution  as  a 
Function  of  Geomagnetic  Latitude 


Figure  3.  Data  Distribution  as  a 
Function  of  Local  Time 
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Figure  4.  Data  Distribution  as  a 
Function  of  Geomagnetic  Activity 
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Figure  5.  Data  Distribution  as  a 
Function  of  Solar  Flux 
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Figure  6.  Data  Distribution  as  a 
Function  of  Altitude 


the  number  of  measurements  durinq  which  the  daily  solar 
flux  value  (with  one  day  laq)  was  within  specific  5  solar 
flux  unit  intervals  is  given.  Figure  6  illustrates  the 
number  of  points  in  each  5  km  altitude  file. 

An  example  of  the  data  base  is  qiven  by  showinq  results 
obtained  at  180  km.  Figure  7  shows  the  density  data, 
corresDondinq  solar  flux  values,  ratio  of  measured  density 
to  the  J71  model  and  the  Kp  index  presented  as  a  function 
of  calendar  year  and  GMT. 


Model  Development 


2 • 1  J71,  J77  and  MSIS  Comparisons 

Histograms  of  the  measured  density  values  were  pro¬ 
vided  to  AFGL  to  permit  analysis  in  relation  to  three 
atmospheric  models  -  J71,  J77  and  MSIS.  Graphical  repre¬ 
sentations  of  frequency  distributions  of  percent  differ¬ 
ences  from  the  models  were  developed.  Figures  8-10 
display  these  histograms  for  the  entire  data  base  for  each 
atmospheric  model.  Figures  11  -  13  display  these  histo¬ 
grams  for  200  km  data  only.  Frequency  distributions  of 
the  data  were  described  in  terms  of  the  mean  value  and  the 
second,  third  and  fourth  moments  about  the  mean. 


2 . 2  Empirical  Model 

A  multiple  linear  regression  analysis  program  was 
utilized  (3)  to  support  the  AFGL  development  of  an  empir¬ 
ical  model  of  atmospheric  density  variations  in  terms  of 
latitude,  longitude,  local  time,  geomagnetic  activity. 


solar  flux,  season  and  semi-annual  effects.  This  program 


Figure  7.  Data  Base  at  180  km:  Solar  Flux,  Density 
Ratio  to  .171  Model  and  Kn  as  a  Function  of  Time 
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Figure  8.  Frequency  Distribution  of 
Percent  Differences  of  All  Data  from 
J71  Model 


Figure  9.  Frequency  Distribution  of 
Percent  Differences  of  all  Data  fn 
MSIS  Model 
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Figure  10.  Frequency  Distribution  of 
Percent  Differences  of  all  Data  from 
J77  Model 
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Figure  11.  Frequency  Distribution  of  Figure  12.  Frequency  Distribution  of 

Percent  Differences  of  200  km  Data  Percent  Differences  of  200  km  Data 

from  J71  Model  from  MSIS  Model 


-•0  -SO  -40  -zo  o  to  to  so  so 


rtsctNT  oirrenrsce  f*on  hbosl  cj??  > 

Figure  13.  Frequency  Distribution  of 
Percent  Differences  of  200  km  Data 
from  J77  Model 
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was  expanded  to  have  the  capability  to  include  a  much 
larger  data  base  -  specifically,  the  composite  AE-C,  AE-D, 
AE-E  and  S3-1  data  base  -  and  more  independent  variable 
terms.  We  then  further  expanded  this  proqram  to  have 
the  capability  to  perform  stepwise  linear  regression  anal¬ 
yses  on  the  data  base  to  determine  the  amount  of  variation 
in  the  data  which  is  described  by  each  term,  the  signifi¬ 
cance  of  each  term  and  the  order  which  each  term  enters 
the  model.  Under  the  present  configuration  the  multiple/ 
stepwise  regression  analysis  program:  (a)  has  the  capa¬ 
bility  to  process  a  data  base  at  least  ten  times  larger 
than  the  original  AE-C  data  base;  (b)  contains  over  one 
hundred  (100)  independent  variables  from  which  may  be 
chosen  the  variables  to  be  used  in  the  particular  model 
to  be  run;  and  (c)  allows  up  to  fifty  independent  vari¬ 
ables  per  model. 

An  expanded  version  of  the  empirical  model  being  de¬ 
veloped  in  support  of  AFGL  contains  21  independent  vari¬ 
ables  and  was  run  at  each  5  km  altitude  from  140  to  240  km. 
The  multiple  regression  program  was  first  run  to  calculate 
the  phases  of  the  trigonometric  terms,  and  then  the  step¬ 
wise  regression  was  run  to  determine  the  order  and  signi¬ 
ficance  of  each  term  in  the  overall  model.  Table  2  dis¬ 
plays  statistical  data  related  to  the  regression  analysis 
including  regression  coefficients,  stepwise  correlation 
coefficients,  and  the  significance  and  order  that  each 
term  entered  the  analysis.  Although  calculated  for  each 
5  km  interval,  these  values  are  given  in  20  km  intervals 
for  this  report. 
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Table  2.  Regression  Analysis  Results 
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2 . 3  Results 

Utilizing  the  regression  analysis  we  were  able  to 
determine  and  display  amplitude  and  phases  of  annual  and 
diurnal  variations  (and  their  harmonics) .  Figures  14  - 
15  display  some  of  these  results  as  a  function  of  satel¬ 
lite  altitude. 

Variations  due  to  semi-annual  effects,  Ap  values,  and 
solar  flux  information  are  given  for  200  km  data  in 
Figures  16  -  18.  Each  figure  displays  data  for  latitudes 
of  0  and  60  degrees. 

Because  the  empirical  model  is  based  on  data  at  selec¬ 
ted  altitudes,  it  is  necessary  to  perform  some  form  of 
interpolation  or  fitting  to  obtain  a  continuous  relation 
between  density  and  altitude.  We  have  chosen  to  fix  (i.e., 
hold  constant)  all  the  parameters  in  the  model  except 
altitude.  This  generates  a  density  value  at  each  altitude, 
and  then  these  values  are  fit  with  a  function  of  the  form: 

d(h)  =  ( c o  +  cjh  +  c2h  +  c3h  )  e  ^a°  + 

where  d  is  the  density  at  an  altitude  h. 

Figure  19  displays  the  resultant  model  fit  as  a  function 
of  altitude  for  March  21  (day  80)  and  June  21  (day  172) . 

The  difference  between  the  two  curves  is  a  consequence  of 
the  change  in  day  of  year.  The  scale  height  H  (=  -1/ai) 
changes  from  25.45  km  for  day  80  to  24.95  km  for  day  172 
which  is  consistent  with  the  larger  range  of  density  shown. 
Table  3  displays  the  model  coefficients  for  the  curves 
displayed  in  Figure  19. 
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Figure  14.  Amplitude  and  Phases  of  Annua]  and  Semi-annual  Variations 
of  the  Data  Base 
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Figure  15.  Amplitude  and  Phases  of  Diurnal  and  Terdiurnal  Variations 
of  the  Data  Base 
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Figure  16-  Semi-Annual  Density  Variation 
for  each  model  at  200  km  for  0  degrees 
and  f>0  degrees  latitude 
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Figure  18.  Solar  Flux  Variations 
of  Density  for  each  model  at 
200  km  for  0  degrees  and  60 
degrees  latitude 
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In  order  to  support  studies  of  density  variation  as  a 
function  of  both  latitude  and  day  of  year,  a  set  of  con¬ 
tour  plots  was  generated.  To  accomplish  this  all  parameters 
except  latitude  and  day  of  year  were  fixed.  Figures  20 
and  21  display  density  contours  as  a  function  of  latitude 
and  day  of  year  for  Ap  =  0  and  Ap  =  80,  respectively.  The 
density  values  are  shown  in  units  of  the  ratio  of  the  em¬ 
pirical  model  density  to  the  actual  average  annual  density 
for  that  altitude,  with  a  contour  interval  of  0.025  in  those 
units . 


3.  Other  Contract  Functions 


3 . 1  EUVS  Experiment  Support 

Utilizing  our  knowledge  (gained  in  supporting  the 
MESA  accelerometer  experiment)  of  the  AE  data  system,  the 
AE  Sigma-9  computer  system,  the  AE  Data  Management  Facility 
(DMF )  software,  and  the  general  concepts  and  philosophy 
of  the  AE  Aeronomy  Team  of  scientists  we  were  required  to 
support  the  efforts  of  the  solar  ultraviolet  experiment 
(EUVS)  scientists  at  AFGL.  We  utilized  the  remote  terminal 
dedicated  to  EUVS  to  provide  interactive  on-line  support. 
This  included  maintenance  and  execution  of  interactive  and/ 
or  batch  programs;  copying  various  output  files  to  the 
remote  printers;  promotion  and  demotion  of  telemetry, 
orbit/attitude,  geophysical  unit  and  unified  abstract  files; 
and  general  account  and  file  maintenance  to  insure  smooth 
execution  of  standard  EUVS  processing.  In  addition,  manual 
plotting,  hand  calculations  and  file  maintenance  programs 
were  written  to  support  this  requirement. 

EUVS  geophysical  unit  data  wore  written  to  magnetic 
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Figure  21.  Density  Contours  as  a  Function  of 
Latitude  and  Day  for  Ap  =  80 
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tape  on  the  Sigma-9  computer,  and  these  tapes  were  then 
shipped  to  AFGL  to  be  read  by  the  CDC  6600  computer.  We 
evaluated  various  tape  formats  and  wrote  programs  to  deter¬ 
mine  the  readability  of  these  tapes  for  the  EUVS  scientists. 

3 . 2  AE-C  Re-Entry  Data 

During  the  re-entry  phase  of  the  AE-C  mission  important 
accelerometer  density  data  were  gathered.  The  data  were 
taken  over  the  final  twelve  days  at  a  rate  of  one  orbit 
per  day,  alternating  spinning  and  despun  mode  data.  Seven 
spinning  and  six  despun  orbits  were  logged. 

Of  the  seven  spinning  orbits  two  orbits  contained  less 
than  three  minutes  of  usable  data;  on  two  orbits  the  in¬ 
strument  was  turned  off;  we  were  able  to  process  data  from 
another  orbit;  and  two  orbits  were  successfully  processed 
to  density  values. 

Of  the  six  despun  orbits  one  orbit  contained  less  than 
three  minutes  of  usable  data,  and  the  other  five  orbits 
were  successfully  processed  to  density  values.  Instrument 
bias  values  were  calculated  from  previous  high  altitude 
circular  despun  data. 

These  re-entry  data  will  be  incorporated  into  the  AE/S3 
density  data  base. 

3 . 3  Circular  Orbit  Processing 

The  obiective  of  the  AE  mission  was  to  study  pheno¬ 
mena  in  the  atmosphere  at  altitudes  above  120  km.  This  was 
to  be  accomplished  in  two  orbital  phases:  elliptic  orbit 
and  circular  orbit.  After  the  initial  elliotic  orbit  phase 
was  completed,  AE-C  and  AE-E  were  circularized  for  the 
second  phase  of  their  mission. 

For  circular  orbit  spinning  data  our  processing  tech¬ 
niques  to  calculate  density  (Reference  (2))  were  for  the 


most  part  unchanged.  For  despun  circular  orbits,  since 
instrument  bias  values  are  required  and  the  satellite  alti¬ 
tude  was  below  the  "non-draq  regions"  used  in  the  elliptic 
orbit  phase,  an  alternate  method  was  required  to  determine 
bias  values.  We  utilized  the  techniques  described  in 
Reference  (5)  to  calculate  bias  temperature  coefficients 
for  each  satellite.  Circular  orbit  temperature  data  was 
then  used  to  calculate  instrument  bias  values,  and  filtering 
techniques  were  used  to  determine  density  values.  Circular 
orbit  density  data  are  stored  in  the  AE  unified  abstract 
files . 

4.  CONCLUSIONS 

An  extensive  data  base  has  been  developed  utilizing 
accelerometer  measurements  obtained  with  four  low  altitude 
satellites.  Utilizing  this  data  base  an  empirical  model 
of  neutral  atmospheric  density  has  been  developed  by  mul¬ 
tiple/stepwise  linear  regression  techniques.  It  is  anti¬ 
cipated  that  these  techniques  and  results  will  aid  in  the 
development  of  improved  models  of  the  lower  thermosphere. 
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